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I . INTRODUCTION 


To meet projected aviation requirements NASA's Langley Research 
Center is conducting a research and development effort, the Terminal 
Configured Vehicle (TCV) program, to develop advanced technology for 
improved terminal area operational capability and safety of transport 
aircraft. The TCV program is aimed at developing capabilities for 
increased terminal area capacity, safe and accurate flight in adverse 
weather conditions including shear winds, noise reduction, the avoidance 
of wake vorticies and reduced fuel consumption. Advances in digital 
flight computers and modern control theory, coupled with accurate 
guidance information such as that provided by the Microwave Landing 
System can be effectively used to achieve some of these goals. The 
work described in this report has evolved within the TCV program frame- 
work [l], [2]. 

The design of the Digital Integrated Automatic Landing System 
(DIALS) was completed in 1980. The automatic landing system was imple- 
mented on the TCV research aircraft, a Boeing 737-100. DIALS was 
flight tested with successful test results by NASA's Langley Research 
Center at Wallops Island Center, thus demonstrating the application 
of modern control theory to a complex design problem. The DIALS 
flight tests were completed in December 1981. DIALS is the first digi- 
tal automatic landing system designed with a modern control structure 
and methodology which has been successfully flight tested, to the au- 
thor's knowledge. This report describes the design and development of 
the Digital Integrated Automatic Landing System prior to flight testing. 



The implementation and flight results will be given in detail in a 
subsequent report. 

DIALS uses the Microwave Landing System (MLS) which is a guidance 
system providing high accuracy position information in the form of 
azimuth, elevation, and range measurements. As the Microwave Landing 
System is less sensitive to weather conditions than conventional 
systems, automatic landing systems using the MLS can be used to reduce 
the congestion in terminal areas due to adverse weather conditions. 
Furthermore, the volumetric coverage provided by the MLS enables the 
use of curved flight paths and steep glideslopes in the final approach 
and landing phase of the flight. The guidance system consists of the 
DME providing range information, an azimuth antenna generally co- 
located with the DME antenna providing the aircraft’s azimuth angle 
relative to the runway centerline up to ±60°, and an elevation antenna 
located at the glidepath intercept point but offset to the side of the 
runway providing the aircraft's elevation angle up to 20°. An onboard 
MLS receiver provides high accuracy position information that can be 
used for steep approaches and curved flight paths in the terminal area. 

This report describes the design and development of a Digital 
Integrated Automatic Landing System (DIALS) for the TCV research air- 
craft, a B-737-100 aircraft through the use of modern control method- 
ologies and structures. The system uses MLS position information, as 
well as onboard sensor measurements; the system was developed using 
modern digital control methodologies [3]. D]> [5], [6]. The phases 
of the final approach and landing considered are localizer and glide- 
slope capture and track, crab/decrab, and flare. The system captures, 
tracks, and flares from a steep glideslope selected by the pilot prior 
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to engaging the capture mode. Thus, the control system modes designed 
are: 

1) Localizer capture 

2) Steep Glideslope capture 

3) Localizer track 

A) Steep Glideslope track 

5) Decrab 

6) Flare 

Two important features of the system are: 

• the simultaneous capture of the localizer and glideslope when 
necessary, and 

• the selectable steep glideslope (2.5° - 5.5° for the B-737) 

Typically, the aircraft is positioned, manually or automatically, below 
the selected glideslope so as to intercept the localizer at a selected 
airspeed, as shown in Fig. 1. When the localizer or glideslope capture 
criteria are satisfied, the corresponding capture mode is engaged; so 
that depending on the initial aircraft position and attitude, the 
localizer and glideslope can be captured simultaneously or sequentially. 
As the aircraft reaches the localizer or glideslope, the corresponding 
track (or hold) modes are engaged. When a cross-wind component is 
present, the control system crabs the aircraft into the wind, followed 
by a decrab maneuver when the decrab altitude is reached. A flare 
path which depends on the glideslope selected for a given approach is 
generated on-line, and the aircraft is controlled about this path until 
touchdown. 
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The main considerations in the development of DIALS include : 


• low overshoots of the localizer and steep glideslope under 
adverse wind conditions 

• quick settling on the trajectory 

• overall smoothness of the flight during the final approach 

• accurate tracking in adverse weather conditions including 
gusts and shear winds 

The ability to follow various steep glideslopes (selectable until 
glideslope capture) provides a flexibility which can increase the 
efficiency of terminal area operations, reduce the noise perceived on 
the ground, and can be used for vortex avoidance when following a large 
aircraft, while reducing fuel consumption during a steep final approach. 
Capturing the glideslope and localizer simultaneously and with quick 
settling times allows close-in captures, while low overshoots of the 
localizer under adverse wind conditions enhance the independence of 
close parallel runway operations. Finally, the low degradation of the 
MLS information accuracy in adverse weather conditions enhances perfor- 
mance capabilities under low visibility conditions. 

The overall objective in the development of DIALS has been the 
direct-digital-design (using modern control methods) of an automatic 
landing system which controls the aircraft from the localizer and glide- 
slope angles between 2.5° - 5.5°, nominal airspeeds between 115 - 135 
knots, aircraft weights between 70,000 - 90,000 lbs, c.g. locations 
between .2 - .3 under adverse wind conditions and low visibility. As 
the control objectives during different phases of the final approach 
vary, the most important desirable characteristics of the control law 
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in the various portions of the flight are different. For example, while 
low overshoot characteristics are important during capture, vertical path 
accuracy becomes more significant in flare. A more detailed description 
of the desirable characteristics in the various control modes is given in 
Section IV. To achieve these varying objectives, it appears that some 
changes in the control law, from one phase of the flight to the next, are 
are necessary. From a strict optimal control point of view, it is 
necessary to change all the control gains for each phase of flight and 
aircraft condition. To avoid many changes, the following approach was 
adopted : 

1) obtain an optimal control law for a nominal condition as a 
starting point, 

2) using realistic non-linear simulations, modify and update this 
law to account for the various non-linearities in the actuator 
servomechanisms, hydraulic systems, engine dynamics and the 
aerodynamic response, as well as to enhance desirable charac- 
teristics which may not adequately be reflected in a quadratic 
cost function, 

3) identify a small number of modifications which can be made at 
each phase of flight to achieve the differing objectives of 
each mode. 

The employment of modern digital design techniques is well-suited 
to the discrete nature of the MLS position information and the use of 
digital flight computers. From the point of view of using lower sampling 
rates, it is also preferable to the approach of designing an analog 
system whose response is then approximated by a digital system. Thus, 
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the control law was developed by imbedding the problem of following a 
specified flight path into the quadratic regulator with disturbances. 

Section II describes the state space formulation of the aircraft 
dynamics and wind models used in the development. Section III formulates 
the problem as a sampled data regulator with disturbances. Section IV 
describes the various modes of the system corresponding to the various 
phases of the final approach and flare. Section V describes the sensors 
used in the filter formulations. Section VI describes the results 
obtained from non-linear simulations of the aircraft and the automatic 
control system under various wind conditions. 
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II. AIRCRAFT DYNAMICS AND WIND MODELS 


The plant model used to develop DIALS was obtained from the 
perturbation equations of the aircraft dynamics about the desired or 
selected glideslope and the localizer at the desired airspeed. The 
wind model and several other state variables were then coupled with 
the perturbation equations to obtain the complete model. The latter 
model was used to generate the control law by applying optimal sampled- 
data control techniques. 

During most of the phases of flight considered, the aircraft main- 
tains a level wings attitude with small deviations of bank angle. Thus, 
the flight condition corresponding to the tracking of the desired glide- 
slope and the localizer was selected as the reference condition about 
which the perturbation equations were obtained. In this condition, the 
longitudinal and lateral perturbation equations are decoupled; i.e., 
small changes in the lateral variables affect the longitudinal variables 
only as second order effects; conversely, the effect of the longitudinal 
variables on the lateral motion is also of second order. Thus, the 
control and modeling of the longitudinal and lateral dynamics were con- 
sidered separately. 

The notation for the coordinate axes used in this study can be 
established by considering three sets of coordinate axes: the earth 

fixed axes, the body axes, and the stability axes. The earth-fixed 
coordinate frame (x^, y^, z g ) has its origin fixed at the glidepath inter- 
cept point (GPIP) on the runway. The x^ axis is along the runway center- 
line, the direction in which the aircraft lands being chosen positive 
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along x^. The axis is the local vertical, positive downwards; y g 

is perpendicular to both x g and Z £ , with its positive end directed so 

as to make the coordinate frame right-handed. The earth is assumed 

to be stationary with respect to inertial space; so that the earth- 

fixed axes form an inertial frame. 

The body axes (x, , y, , z, ) and the stability axes (x , y , z ) are 
bbb sss 

fixed to the body of the aircraft; i.e., they are body-fixed axes. The 

origin of both axes is fixed at the aircraft center of mass. The x, 

b 

axis is along the fuselage reference line of the aircraft, positive 

towards the nose, the y^ axis is positive towards the tip of the right 

wing, the z^ axis is perpendicular to both x^ and y^ and is positive 

downwards (when the aircraft pitch angle is zero). The stability axes 

(x , y , z ) are obtained from the body axes by a rotation of a , the 
sss o 

steady state angle of attack, about the y^ axis. The three sets of 
coordinate frames are shown in Figure 2. 

A. LONGITUDINAL EQUATIONS 


The general equations of motion for rigid aircraft can be linearized 
about a steady flight condition as described in [ 7 ], [ 8 J. The nominal 
flight condition used here corresponds to flight among the selected glide- 
slope and localizer, at a constant airspeed with the flaps at 40 degrees 
and the gear down in the landing configuration. The longitudinal pertur- 
bation equations for describing the aircraft's motion in the vertical 
plane can be expressed in the stability axes as: 

mu = - rag cosY q 6 + f & + f T (1) 

x x 
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( 2 ) 


m(w - U q) = - mg smY 0 + f . + f 

o o AT 

z z 


I 

yy 


where 


U 

o 



I 

yy 


m 

u 

w 


q = m A + n, T 


(3) 


- steady inertial speed in the x g direction 

- glideslope angle 

- perturbation in pitch angle 

- pitch rate 

- perturbation in net aerodynamic force along the x g direction 

- perturbation in net aerodynamic force along the z direction 

s 

- perturbation in thrust along the x direction 

s 

- perturbation in thrust along the z direction 

s 

- perturbation in pitching moment due to aerodynamic forces 

- perturbation in pitching moment due to thrust 

- moment of inertia about the y axis 

y s 

- aircraft mass 

- perturbation in inertial speed along the x direction 

s 

- inertial speed along the z^ direction 


In equation (1), the term f represents the total algebraic change 

x 

in the value of the aerodynamic force along the x g axis due to changes 

from steady state values in the values of the aerodynamic and control 

variables; the terms f^ , f^, , f^ , m^ and m^ are defined similarily as 

z z x 

the changes in the appropriate forces or moments from their steady values 
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on the glideslope. These terms can be expressed in terms of the aircraft 
stability derivatives, the moments of inertia and the perturbations in the 
aerodynamic and control variables. Substituting these expressions for the 
force and moments into equations (1) - (3) , the aircraft equations of motion 
can be expressed as given below. 


mu = - mg cosY q 0 + q Q S |(" C D + + C T + 2C T ) i± T 

* u' O XIJ ' xo 

+ (C T - CL ) a - C_ * 6e - C x 6s + (L, <5T . > 

Lo Da — Doe Dos T x oT ) 


(4) 


m(w - U q) = - mg siny 0 + q S < - (CL , 4- 2 CL ) u T 
o o o ( Lu_ Lo — 

- (C T + CL. ) a - CL . a - CL q l 
Ly ^ = 5 0 Sc{ + 2 C Mo ) u' + (C^ + a + a 

+ Sla 4 + C M6e 6e + C MSs & + C M«T 5T } ’ 

where 


(5) 


( 6 ) 


u + u 
» W 


u = 


U 


(7) 


a = a + a , 
— w 


^ = q + q w 


( 8 ) 

(9) 


u w > a w > and q^ are the components due to wind, q^ is the steady value of 
the dynamic pressure at the selected airspeed, S is the effective wing area, 
6e, 6s, and 6T are the perturbations of elevator, stabilizer and thrust, 
respectively. 
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Equations (4) - (6) describe the linear and angular velocities of the 
aircraft in the vertical plane. The position of the aircraft can be obtained 
by integrating the inertial velocity components over time. Thus, 


x 

= tT = l es (1,1) (1 + u,) + l es (1>2)S + L ES (1 ’ 3) “ ’ (10) 

o 

2 

*4 = ^ = L es (3,1) (1 + u’) + L es (3,2)6 + L es (3,3)a , (11) 

O 

where L ££ (i,j) is the element in i* 1 * 1 row and j column of the matrix, 
Lgg, which represents the transformation from the stability axes to the 
earth-fixed axes. Note that a and $ are inertial quantities and corres- 
pond to normalized velocity components in the stability axes; under no 
wind conditions these would be the same as the aerodynamic angle of attack 
and sideslip. The position equations can be rewritten in the following 
form. 


x T = -sinY 
e o 

e 

+ cosy u* 
o 

+ sinY 

o 

a + n x , 

(12) 

z ’ = -cosy 
e o 

0 

- sinY u f 
o 

+ cosy 

o 

a + n z , 

(13) 


where 


n x = l es (1»1) (1 + u') - cosY q u 1 + siny Q 9 + L £ g (1,2)3 

+ (L ES (1,3) - sinY Q )a 

n z = l es (3,1) (1 + u r ) + sinY Q u f + cosY q 0 + L £S (3,2)3 

+ (L £S (3,3) - cosY Q )a 


(14) 


(15) 
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In this form, the differential equations are linear with a forcing 
function that contains the non-linear part, which are second order terms 
with respect to the steady flight conditions considered. With this 
approach we can use linear theory in the development of the filter and 
control law without neglecting the non-linear terms completely. To 
obtain the second order terms for the remaining variables, consider the 
equations 


0 = cose}) q - sincj) r , 

a = u + qw - rv , 
x 

s 

a = w - q(U + u) + pv 
z o 

s 

Rearranging these equations, we obtain 


(16) 

(17) 

(18) 


0 = q + n ( 


(19) 


u = 


ir + V 

o 


a = 


w 


U 


+ q + n 


a ’ 


( 20 ) 


( 21 ) 


where 


r| Q = (cos(j) - 1) q - sin<£ r , (22) 

n 1 = - q a + r 3 , (23) 

r| = q u' - p 3 (24) 

a 

It may be noted that the perturbation equations (4) - (6) neglect 

the second order terms represented by ri , n , and n • The second order 

u u cx 
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terms introduce some coupling between the longitudinal and lateral equations. 
To model the actuator dynamics and obtain a control rate structure 


6T = -.5 6T + .298 6th , 

(25) 

6th = u^ , 

(26) 

6s = u 3 

(27) 


where ST is the thrust perturbation in units of one thousand pounds per 
unit of ST, Sth is the throttle perturbation in degrees and Ss is the 
stabilizer perturbation in radians. As the lags in the elevator action 
are small, the elevator time constant was neglected. The aircraft’s 
longitudinal equations of motion, the position equations and the actuator 
equations can be combined and after some manipulation can be expressed in 
state variable form. 


w, 


= (0 

u' a q x’ z' 6 t 6th 6s) 

> 

(28) 

= (6e 

6s 6th) , 


(29) 

= (u* 
v w 

a q ) 
w w 


(30) 

" (n e 

n T nOrir|0 0 0) 

u’ a x z 

J 

(31) 


*p = A.x, + B u ? + DpW p + rip , 

A/ A/ Aj Aj Aj A j Aj Aj 


(32) 


where A^, B^, and are matrices of appropriate size corresponding to the 
coefficients in the original equations. Expressions for the elements of 
these matrices can be obtained in terms of the stability derivatives of 
the aircraft. 
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To complete the model of the longitudinal equations of motion, it is 
necessary to model the wind velocities which affect the motion of the air- 
craft as seen in (32) . 

The longitudinal wind model contains the components of steady wind 

velocities, turbulence and shear winds in the longitudinal axes. The 

turbulence model uses the Dryden spectra ^4j] for the various components 

varying with altitude. The turbulence model has three components: u’ in 

S 

the x. direction, a in the z, direction, and q which models the effect 
b g b g 

of turbulence on the pitch rate of the aircraft. These components are 
modeled using the following spectra. 


S (ft) 
u 


20 2 L 
u u 

1 + (L ft) 
u 


S (ft) 
a 


G 2 L 1 + 3 (L ft) 2 
w w w 


2 2 2 ’ 
V 1 + (L ft) 

a w 


S (ft) 

q 


2 2 
JT V 

a 


1 + (— ft) 2 a 

TT 


s (ft) 


(33) 


(34) 


(35) 


where b is the wing span, L and L are the scales of turbulence, V is 

& r u w a 

the airspeed, and ft is the spatial frequency related to the temporal 
frequency co by 


ft = w/V . (36) 

a 

The u' component is independent of a and q ; however, a and q are corre- 
g g g § S 

lated with their cross spectral density being 
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03 


( 37 ) 


S ( 03 ) = 
qa 


jw 


1 + j 


4b 

ttV 


S ( 03 ) 

a 


The above spectra can be factored using spectral factorization methods to 
obtain a linear system driven by white noise which generates an output having 
the above spectral characteristics [[ 7 ], [ 3 J, [ 13 j. Thus, the following 


transfer functions are obtained to generate u', a and q , 

■ 8 8 6 


G u (s) = 


1 s 

a 


(38) 


G a (s) = 


, , , w 

1 + 3 — s 

a 


L L „ 9 

1 + 2 — s + (— ) s 

a a 


(39) 


G q (s) “ , . 4b 


(40 


1 + 


TTV 


where a is the input to the system G (s) to obtain q with the specified 

$> q s 

spectrum and cross-spectral density. Figure 3 shows block diagrams of 
the systems generating the turbulence components. 

The steady and shear wind in the longitudinal direction was modeled 
by 


u = u- 
s sh 


sh 


= s 


£3 


(41) 


a 

s 



(42) 


Thus, to simulate a specified shear profile for u^, with appropriate 

initial conditions, e.g., to obtain a linear profile u^ changing at a rate 

of u\ , the initial conditions for u f , is set to u 1 , and is set equal 

sho sh sho £3 
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to zero; alternately, an impulse in E,^ will also achieve the same profile. 

The transfer functions obtained for the wind model can equivalently 
be expressed as differential equations in state variable form as shown in 
(43), 


W a A wJ, W £ + B w£ 


w & = C w5, V SL 


where w„ is given by equation (30). Thus, the coupled equations (32) and 

X/ 

(43) model the longitudinal motion of the aircraft under various wind 
conditions. 


B. LATERAL EQUATIONS 

The lateral equations of motion were obtained using the same approach 
as the longitudinal equations. The flight condition about which the general 
equations of motion were linearized corresponds to flight along the local- 
izer, at the selected glideslope and airspeed, with the flaps at 40 degrees 
and the gear down in the landing configuration. For this condition, the 
perturbation equations can be written as 


m(v + U r) = mg cosy 4> + f + F , 

° ° -y 1 * * * y 


I p-I r = JL + £ » 

xx r xz AT 


I r - I p = b.+n_ 

zz xz A T 


(44) 

(45) 

(46) 


Equations (44) - (46) are the perturbation equations which describe the 
lateral dynamics of the aircraft, where 

U - steady-state inertial speed along the stability x-axis, i.e., 

along x , 
s 
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v - perturbation in the inertial speed along the stability y-axis, 
i.e. , along y g , 

p - perturbation in the roll rate about x^, 
r - perturbation in the yaw rate about z^, 

Y q - steady-state flight path angle, 
cj) - perturbation in the roll angle, 
f - perturbation in the net aerodynamic force along y , 

Ay s 

f - perturbation in the thrust force along y , 

i y . * S 

£^ - perturbation in the rolling moment due to aerodynamic forces, 

n^ - perturbation in the yawing moment due to aerodynamic forces, 

- perturbation in the rolling moment due to thrust, 

n - perturbation in the yawing moment due to thrust, 

and I , I and I are the moments of inertia of the aircraft in the 
xx zz xz 

stability axes. 

The terms f , f , £ , £ , n. and n can be expressed in terms of the 

Ay -L y A _L A _L 

stability derivatives of the aircraft evaluated at the steady state values 
of the aerodynamic variables and the control surface settings, in linear 
form [?].[ 8], [9]. Substituting these expressions into equations (44) - 
(46) and rearranging terms, we obtain linear differential equations in the 
sideslip angle, the roll rate and yaw rate. Writing the derivatives of 
the roll and yaw angles in terms of roll and yaw rates, we obtain the 
following set of differential equations. 



(47) 

(48) 

(49) 
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p = a 43 6 + a 44 P + a 45 r + t> 41 6A + b 42 6R + b 43 Ssp 


+ h , 1 3 + h. 0 

41 w 42 


P + 
w 


43 


r 

w 


r = a co B + a_, p + a cc r + b__ 5A + b co 6R + b co 6sp 

53 54 55 51 52 53 


+ h 


51 


B + h_ 0 p + 
w 52 w 


53 


r 

w 


( 50 ) 


(51) 


where B s p and r are the sideslip angle, roll rate and yaw rate, respec- 
tively, B > P and r are the sideslip angle, the roll rate and yaw rate 

WWW 

due to wind velocities only, 6A, 6R and 6sp are the perturbations of the 
ailerons, rudder and spoilers, respectively. The coefficients in the 

above equations depend on the aircraft stability derivatives and are given 
in the Appendix. Thus, a set of linear differential equations describing 
the lateral velocities and attitude of the aircraft can be obtained. 

The position of the aircraft relative to runway centerline is expressed 
by the perpendicular distance of the aircraft center of mass to the runway 
centerline. This distance normalized by the aircraft’s steady-state speed 
will be used as a state variable in addition to the equations already ob- 
tained. Then, the lateral distance y (in feet) of the aircraft can be 
expressed by 


y = U o [L es ( 2,1)(1 + u') + L es ( 2,2)6 + L ES (2,3)a] , (52) 

where u' is the normalized inertial perturbation in the speed along the 
x g direction and a is the perturbation in the inertial angle of attack, and 

L__(2,l) = cosot cos (9 + 0) sinijj + coscj) sina sin (9 + 0) sin\{> 

ES o o o o 

- sina Q sin(f) cosip (53) 
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( 54 ) 


L„„(2,2) = sin(J) sin(0 + 0) siniJJ + cos<|> cost|> 

Eib O 


L wt? (2»3) = -sina cos (0 +0) simii + cosa cosd) sin(0 + 0) sinib 

lib o o o o 


- cosa^ sin<|> cosip 


(55) 


Rearranging the terms in equation (52) 


y ? = 3 + cos(a - 0 )ip - sina d> + ri 
y o o o y 


(56) 


where T] and y 1 are given by 


n y = (L eS (2>2) - 1)3 + L es ( 2,1)(1 + u f ) - cos(a^ - djip 


o o 


+ L (2,3)a + sina (J> , 

JtjiJ o 


(57) 


y ' = y/u 


(58) 


Note that equation (56) contains no approximation when a in (57) is 
interpreted as the normalized inertial velocity component in the z g 
direction, but is simply a rearranged form of (52) with the non-linear 
terms grouped into a single term. To obtain the second order terms for 
the remaining variables, note that 


d> = sec© (cosy p + sinv r) + r\, , 
o ' o ' o <p 


ip = sec0 (sina p + .cosa r) + n, , 
r o o o ip 


(59) 

(60) 


6 ■ - r + i e * 


(61) 


where 


= (cos<f> tan(0^ + 0) - tan0^) (cosa^ r + sina^ p) , (62) 
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( 63 ) 


r\, - (cos<p sec (6 + 0) - seed ) (cosa r + sina p) , 

(p o o o o 

rig = pa - ru , (64) 

when q is assumed to be small. 

The servo systems, hydraulic and mechanical actuator systems on this 
aircraft have clearly noticeable non-linear effects. These effects include 
usual non-linearities such as hysteresis, rate and position limits, as 
well as d.c. gain variations as a function of actuator position, dynamic 
pressure, etc., and non-linear spoiler feedforward and feedback systems. 

These non-linear effects were not included in the ’’control design model” 
described here. However, these non-linearities were included in the simu- 
lation model, and design changes were made to accommodate these non-linear 
effects during the control development using the non-linear simulation and 
later during the flight tests. 

The controls which affect the lateral motion of the aircraft are the 
aileron, rudder and spoiler surface settings as can be seen from equations 
(47) - (51), where the spoiler action is used mainly to aid the effect of 
the ailerons during turns. Thus, the spoiler setting is programmed according 
to the aileron setting. This is approximated and modeled here as 


6s = C 6 A 

p spa 


C = 1.73 . 

spa 


(65) 


A rate command structure was used for the rudder control; hence, the 
rudder position is considered a state variable which is obtained by inte- 
grating the rudder rate command. 



( 66 ) 
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where is considered to be the rudder rate control. 

If the relation between the spoiler and aileron given in (65) is 
substituted in (49), (50) and (51), then the spoiler terms are eliminated 
from the equations. Now, forming a state vector x such that 

x T = (cj) ip 8 p r y’ 6R) (67) 

and a control vector u such that 

u T = (6A 6R) , (68) 

equations (47) - (51) , (56) , (65) and (66) can be combined into a state 
variable model of the lateral motion of the aircraft of the form 

x = Ax + Bu + Dw + r| , (69) 

T -T 

where w = (8 p r ) and n = (0. n. n Q 0 0 n 0) 

www (p ip 8 y 

The lateral motion of the aircraft is described by the state variable 

model given in equation (69); this model describes the response of the 

aircraft when a control is applied or when the wind velocities such as 

gusts or steady winds are non-zero. The effects of the wind velocities are 

introduced through the vector w. The components of this vector are 8 or 

w 

the wind velocity along the y g direction normalized by the airspeed of the 

aircraft, p or the rotation of the air around the aircraft about the x 
w s 

axis, and r or the rotation of the air around the aircraft about the z 
w s 

axis, respectively. The roll rate p^ and yaw rate r^ components of the 
wind vector w consist only of the effects of wind gusts, thus having an 
average value of zero; i.e., these components do not have a steady state 
effect but introduce turbulence effects into the equations. On the other 
hand, the 8^ or the normalized lateralwind velocity contains terms for 
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both wind gusts and steady winds; thus, it is modeled as 


3 = 6+3 

W g s 


( 70 ) 


where (3 is the gust or turbulence term, and 3 is the steady wind term. 

g S' 

The gust terms are of a random nature and can be modeled using the well- 
known Dryden spectrum m,[8]. This method consists of using spectral 
factorization methods to obtain a dynamical system which generates a random 
process having the specified power spectral density when driven by a white 
noise process [10], [11], [12 ] , and [13]. 

The Dryden spectra describe the statistical behaviour of wind gust 
velocities along the aircraft body coordinates by specifying their power 
spectral densities in terms of the spatial frequency ft [14], The spectra 
for the gust components of interest are given below. 


V n) - < 


L 1 + 3 (L fi) 2 

V V 

7 TV 2 [1 + (L tt) 2 ] 2 
a v 


S (ft) 
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a 


7TL Wn 1/3 
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4bJ2 2 

7T 


S (ft) 
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1 + ( 



3b ft 2 

TT ' 


Sg(fi) 


(71) 


(72) 


(73) 


where b is the wing span, L and L are the scales of turbulence, V_ is the 

v w a 

2 2 

airspeed, a is the variance of the lateral gust and a is the variance of 
v w 

the vertical gust. The change from the spatial frequency ft to the temporal 
frequency CO can be made by 



( 74 ) 
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Substituting equation (74) into the expressions for the power spectral 
densities of the wind gusts, we obtain the spectra in terms of the temporal 
frequency; then using spectral factorization techniques the following trans- 
fer functions can be obtained 


Gg(s) 


-a 

v 



3 L s 

1 + -vr~ 


(1 + 


ky s ) 


V s 


(75) 


G (s) 
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1 

TT L _ JL 
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i I 2 

, 8 (-4^3 

2 

0 

W 

L V 

i,4 b s 



w a 

1 + 




TT V a 



(76) 


G r (s) 


-s 


1 + 


3b s 

TT 


V s) 


(77) 


It should be noted that even though p^ is independent of 3 w and r^, 

the latter two are not independent of each other. Thus, if a white noise 

process is input to the transfer function G^(s), the output would have the 

desired power spectral density, but may not have the desired cross-correlation 

with 3 • Hence, equation (77) must be interpreted as 3 being the input to 
S § 

the first term in the above equation in order that the proper cross-cor- 
relation be obtained. 

The wind gust terms can thus be simulated by passing white noise 
through the systems with the transfer functions given in equations (75), 

(76) and (77). The lateral wind, however, has a steady or average value 
which is not necessarily negligible. Thus, consider that a steady wind is 
present; in the earth-fixed coordinate system, the wind velocity has a 
component in the direction of runway centerline W^, and a component 
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perpendicular to the runway centerline say W^; it is assumed that there is 
no steady wind in the vertical direction although gusts may be present. 
Hence, if Lg E is the transformation matrix from earth-fixed to stability 
coordinates, then the steady component, $ s> of the normalized lateral wind 
velocity is given by 

■ L SE (2 ’ 2) W 5 + W 2 ’ 15 w x > < 78 > 

Lg E (2,l) = sin(0 Q + 0) sin(f> cosi|; - cos(f> sinijj , (79) 

Lg E (2,2) = sin(0 Q + 0) sin<j> sinip + cos<J> cos ip . (80) 

To include wind shear into the model, the steady lateral wind velocity 
can be described as 


^5 “ W 6 + W 3 » (81) 
W 6 = w 4 , (82) 

where and are gaussian white noise processes independent of each 

other and of 3 , P and r . Now, the transfer function for the gusts 
g w w & 

described in (75) , (76) and (77) can be combined into a state variable 
model of fourth order. Adding (81) and (82) to this model, we obtain a 
sixth order model of the form 


W = A W + B 0) 
w w 


where 3 W can be expressed as 


3 w = 6 g + 3 s ^ W 1 + W 5 + *i 


h = (L SE (2 ’ 2) " 1)W 5 +L SE (2 ’ 1)W x 


(83) 


(84) 

(85) 
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Defining the remaining elements of the vector £ to be zero, w can be written 


as 


w = C W + £ 
w 

This expression can now be subsituted into (69) to obtain 


( 86 ) 


x = Ax 4- Bu + DC W + n (87) 

W 

n = d£ + n (88) 

The model for the lateral dynamics, including the aircraft aerodynamics, 
the actuators and wind conditions, is given by (87) and (83). 
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III. PROBLEM FORMULATION 


The mathematical models obtained in the last section can be used to 
formulate a discrete stochastic optimal control problem for the design of 
a digital automatic landing system. The solution to this optimal control 
problem was used as an initial design in the development of DIALS. 

Although some modifications to the design were made to account for non- 
linearities and other unmodeled effects, the basic structure of the optimal 
control law was left unchanged. 

The design problem was formulated as a discrete stochastic quadratic 
regulator with random disturbances OIL 5 J, £ 6j. The longitudinal and 
lateral control laws were obtained separately, but using the same basic 
approach. As seen from (32), (43), and (87), (83), the state variable 
form of both the longitudinal and lateral models is the same; so that 
either model can be expressed as 

x = Ax + Bu + DC W + n , (89) 

W 

W = A W + B w . (90) 

w w 

To obtain a constant gain control law for the various phases of the 
final approach and landing, the desired trajectory was modeled in the form 


z = A z + £ , (91) 

z ^z 

- s 2 • 

where is assumed to be a white noise process. The desired trajectory, 
z(t), is thus modeled as a random process whose statistical properties are 
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determined by (91), (92). The assumption of a random trajectory is a con- 
ceptually appealing one, as it stresses the fact that the trajectory may 
be altered at any time; the future values of the trajectory are thus not 
certain, but can be predicted using the model of the trajectory. In gen- 
eral, for an arbitrary deterministic trajectory, the optimal control 
depends on future values of the trajectory £ 15 ], except when the trajectory 
can be expressed as the output of a homogeneous dynamic system. 

The error in the actual trajectory, x(t), can now be defined as 

e(t) = x ( t ) - z(t) . (93) 


The models of the aircraft dynamics (89), (90) and the desired trajectory 
(91), (92) can be combined, and expressed in terms of the trajectory 


error, e(t), in the form 

e = Ae + Bu + DC W + (A - A )z 
w z 

W = A W + B a) 
w w 

z = A z + c 
z z 



where 03, E, and £ are assumed to be 
z ri 

Rearranging the terms in (94) - 
e = Ae + Bu + Dd , 
d = A d d + £ , 


- C z + n , (94) 

(95) 

(96) 

(97) 

(98) 

gaussian white noise processes. 

(98), 

(99) 

( 100 ) 
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where 


,T /tt T T _T T.T 

d = (W z c z n ) 


(101) 


and E, is a white noise vector. In (99), (100), the disturbance form of 
the problem is apparent, where the disturbance, d, is uncontrollable by u, 
and is the output of a marginally stable system. 

A cost function quadratic in the error, e(t), of the form 


J c = 2 ~— / E |e T (t) Q e(t) + u T (t) R u(t) |dt 


( 102 ) 


'f o 

was selected for the initial design. To obtain a discrete control law for 
digital implementation, the usual assumption of a constant control between 
sampling instants was made. 


u(t) = u k , kT Q < t < (k+l)T 0 . (103) 

With the assumption of (103), the system equations in (99) and (100), 
and the cost function in (102) can be expressed in terms of the state 
variables at the sampling instants Ol 


e k+i = * e k + ru k + r d d k + 5 ik 


d k+l = d k + 5 2k 


(104) 

(105) 


where e, and d, represent the samples e(kT ) and d(kT ), and C,, cj 0 , 

K K. O O IK ZK 

are white noise sequences. Similarly, the cost function (102) can be 
discretized in the form 
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J = 


2 (N+l) 


" r t 

JoL e k 


5 e k + \ £ u k + 2 4 s d k + 2 % u k 


+ 2 d. s 

k 




(106) 


It is important to note that the discrete cost J and the continuous 
cost J c differ by a constant independent of the control sequence used. Thus, 
to find the control sequence which minimizes the continuous cost J^, it 
suffices to obtain the control which minimizes J. The control which mini- 
mizes the cost J subject to the constraints of (104) and (105), is given 

by [6], 


U k = ~ H lk \ - H 2k d k 


(107) 


where and d^ are the least mean-square estimates of e^ and d^, respectively, 
given past and current measurements, and the gains H^_ and H ?v are given by 


2k 


:-l 


H lk \ G lk 


H 


2k 


\ 


-1 


r 2k 


G ik = r T p lk * + m , c 2k = r T D k + s , 


D k " F ik r a + p 2k \ = S + rT p ik r • 


(108) 

(109) 

( 110 ) 


while P-. and P„, are given by the nonlinear difference equations: 
lk 2k 


r p rn A 1 A A 

P ik-i = ^ P lk*- G ik\ G lk + Q • p in = q ’ 
P 2k-1 " » ' r £ 'k G lk ]T \ + S ’ P 2N - S > 


( 111 ) 

(112) 


Thus, over a finite optimization interval t^ or N, the optimal control 
law is specified by (107) - (112). It should be noted that the gain 
and the cost matrix P^ are the optimal solutions to the LQG problem without 


any disturbance; i.e., d, = 0. As the optimization interval N increases. 
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it is well-known that and P^ ->• P^, under loose conditions []l7], 

where the closed-loop system matrix (p - TH, is stable. If the disturbance 
matrix <f)^ satisfies p(<{>^) - 1» it can be shown that P^^ P 2 as the opti- 
mization interval N increases without bounds, and the cost, J, converges 
to a finite value. It should be noted that, due to the random plant noise 
added at every sample, the averaging of . the cost over the optimization 
interval is ncessary to maintain a finite stochastic cost in any (stochastic) 
LQG problem. The gains and H 9 can then be obtained by solving (108) - 
(112). Using the 
rewritten as 

“k = - H e 
where the control 


se steady state gains, the control law in (107). can be 


- H W - H Z. - H C - H n, 
wk zk £ k Ok 


(107 a) 


gains for each term in d^ are shown explicitly. 
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IV. DEVELOPMENT OF THE CONTROL LAW 


The main objective of the development was to investigate a direct- 
digital-design using modern control techniques that would lead to flight 
tests of the automatic landing system designed on the TCV B-737 Research 
Aircraft . The phases of the final approach and landing considered were 
the localizer and glideslope capture and track, crab/decrab and flare to 
touchdown. The range of flight conditions considered was 2.5° - 5.5° for 
glideslope angle, 115 - 135 knots for airspeed, 70,000 - 90,000 lbs for 
aircraft weight, .2 - .3 for c.g. location, and crosswinds below 15 knots. 

The main considerations in the development of DIALS include the over- 
all smoothness of the flight during the final approach, low overshoots of 
the localizer and glideslope, quick settling on the trajectory and accu- 
rate tracking in adverse weather conditions including gusts and shear winds 
However, the relative importance of these desirable characteristics varies 
with the phase of the final approach considered. For example, low over- 
shoot characteristics under differing wind conditions are particularly 
important during the localizer and glideslope capture maneuvers whereas the 
actual capture path is of secondary importance; however, during the flare 
maneuver, the vertical path error is more significant, and, in particular, 
the touchdown pitch attitude must be sufficiently positive to avoid landing 
on the nose gear. Similarly, when tracking the localizer and glideslope, 
it is important to avoid any offsets and have low sensitivity to wind gusts 
A more detailed description of the desirable characteristics in the various 


The DIALS flight tests were completed in December 1981. The flight data 
and test results will be described in a separate report. 
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control modes is given in the following. To achieve these varying objectives, 
it appears that some changes in the control law, from one phase of the flight 
to the next, are necessary. 

Some of the differing objectives for each phase of the final approach 
can be accommodated by different values of the weighting matrices (Q, R) 
in the cost function. However, in an optimal control design, this would 
necessitate changing all of the gains of the control law at the initiation 
of each mode, along with possible transient problems due to switching gains. 

To avoid possible complexity, the approach taken was to obtain a basic 
design considering the total trajectory, and then vary a small number of 
gains using easy-on’s and to modify the structure by switching integrators 
in or out appropriately, at the initiation of each mode. Thus, the basic 
design model does not contain any integrators; however, when the objective 
requires it, integrators are switched in or out. For example, at the 
initiation of the localizer (or glideslope) track mode, the integral of the 
lateral (or vertical) offset is introduced to obtain a type 1 system without 
steady offsets. Thus, the basic design is obtained as a stable closed-loop 
system to which modifications are switched in or out at each mode according 
to the main objectives of that mode. The specific modifications made at 
each phase are detailed below. 

Non-linearities in the system and unmodeled effects in the aerodynamic 
characteristics, actuator systems, sensors and electronics were simulated 
in as much detail as available to obtain a realistic computer simulation 
of the overall system. These effects were then reduced or eliminated through 
simulation analyses. 
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Thus the approach to the overall design of the digital automatic land- 
ing system (before flight tests) was: 

1) obtain an optimal control law for a nominal condition as a 
starting point, 

2) using realistic non-linear simulations, modify and update this 
law to account for the various non-linearities in the actuator servo- 
mechanisms, hydraulic systems, engine dynamics, and the aerodynamic 
response, as well as to enhance desirable characteristics which may not 
adequately be reflected in a quadratic cost function, 

3) identify a small number of modifications which can be made at 
each phase of flight to achieve the differing objectives of each mode. 

It should also be noted that the cost function (102) penalizes the 
error in the actual trajectory relative to the desired trajectory; however, 
it also penalizes the total control u. As noted in £ 6], such a cost fun- 
ction does not necessarily provide sufficient trim unless the cost matrices 
are appropriately selected or the control cost is modified to penalize the 
deviation from a pre-selected trim [18]. For example, consider the action 
of the elevator and stabilizer. On the TCV B-737 Research Aircraft, the 
effects of the elevator and stabilizer on the motion of the aircraft are 
the same, except that the effectiveness of the larger stabilizer surface is 
twice that of the elevator. It is desirable to make trimming changes using 
the stabilizer, and maintain the elevator near its equilibrium position 
(referenced at 0) to avoid high moments on the elevator hinge over sustained 
periods of time. Of course, the fast motion (due to the actuators) of the 
elevator is used to stabilize the dynamic modes of the aircraft. Since the 
trim values required on different glideslopes, or wind conditions, or c.g. 
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locations, etc. vary with the particular condition encountered, it is de- 
sirable to have the stabilizer position move to trim the aircraft, while 
the elevator is maintained at low hinge moment values except in transients. 
Thus, it is desirable not to penalize stabilizer position, but to penalize 
elevator position to obtain the desired trimming action. Since, it is 
also desirable (and necessary) to have slow stabilizer motion, high stabilizer 
rate should be appropriately penalized to obtain smoother transitions. Thus, 
appropriate selection of the cost matrices with rate command structure can 
provide appropriate trimming with optimal control formulation described in 
the previous section. Similarly, throttle position is not penalized, while 
a non-zero penalty on throttle rate provides smoother and slower autothrottle 
action. In general, control activity can be placed at desirable levels by 
penalizing the control rate. 

To obtain a 3-D control law, the variable corresponding to distance 
along runway centerline (x^,.) was eliminated from the equations of motion 
for the design model. However, the filter development model included (and 
estimated) this variable. 

The specific modification used at the initiation of the various modes 
of the automatic landing system are described in the following sections. 


A. LOCALIZER CAPTURE 


The main function during the localizer capture mode is to perform a 
smooth transition from the aircraft’s initial position to the runway center- 
line, so that the aircraft can quickly settle and stabilize on the localizer. 
Two important characteristics were stressed in the design and development 
of DIALS for this mode. 
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The first characteristic is the reduction of overshoots of the localizer 


and the settling time under adverse weather conditions. Consistently low 
overshoots of the localizer are particularly important in airports with 
close parallel runways. This property can be used to increase the efficiency 
of the terminal area as it allows the two parallel runways to be operated 
independently. A low overshoot also reduces the settling time on the 
localizer, thus allowing close-in captures when necessary. 

The second characteristic is to initiate the capture maneuver by 
rolling the aircraft away from the runway centerline. Pilots usually 
initiate localizer capture in this fashion, and consider rolling towards 
the runway centerline "the wrong direction". Thus, this characteristic 
makes pilot monitoring of the capture easier, as it is the expected behaviour. 
This property can be obtained by initiating the localizer capture when the 
natural tendency of the control law is to roll in the desired direction. 

To achieve this property, consider the component of u in (107a) corresponding 
to the aileron command. The distance, y^, from the localizer at which the 
aileron command crosses zero was selected to initiate the localizer capture 
mode. Assuming level wings and using the initial conditions for the filter 
states, and neglecting small terms, a localizer capture criterion of the 
following form can be obtained 


l*J < U cJ c l + c 2 K 

+ c 3 

6 W J * 

(113) 

C 1 = H e 1 2 /H e 1 6 

C 2 

= H e 1 3 /H e 1 6 

(113 a) 
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where y , ip , (3,, and 3 , are estimated of the distance from runway center- 

Jv K K WK. 

line, the heading relative to the runway, the inertial normalized side 
velocity and the normalized cross-wind velocity, respectively. 

The localizer capture mode is initiated when the inequality holds. 

The constants are chosen so that the initial tendency is to roll away 
from the localizer with the possible exception of conditions when the 
aircraft is already rolled as desired, say due to gusting winds. It 
should be noted that the localizer capture mode is engaged farther away 
from the runway centerline when the desired speed U q is higher, as would 
be expected. To further ensure a smooth transition, an "easy-on" function 
was used on the aileron command signal. 

To avoid large overshoots beyond the runway centerline, appropriate 
closed-loop damping of the control law was obtained by proper selection 
of the weighting matrices Q and R and simulation of the localizer capture. 


B. GLIDESLOPE CAPTURE 


The desirable properties of the glideslope capture mode include low 
overshoot and the tendency to pitch down at the initiation of the capture. 
To have a low overshoot of the glideslope for the various steep glideslope 
angles that may be selected (i.e., 2.5° - 5.5°), it is desirable to ini- 
tiate the capture when the aircraft is well below the glideslope. This 
approach is different than the practice of engaging the glideslope hold 
mode when the glideslope is intercepted (e.g., see [9] ) , where an over- 
shoot, even if small, cannot be avoided. 

Initiation of the capture mode below the glideslope alone does not 
guarantee a low overshoot. Thus, some experimentation in the selection 
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of the weighting matrices coupled with simulation was used to achieve this 
property. In particular, for both overshoot characteristics and smoothness 
of maneuver, the sink rate command was eased on from zero to the glideslope 
sink rate. This implementation provides a simple method of adjusting the 
speed with which the capture maneuver is performed without an appreciable 
overshoot . 

While engaging the capture mode below the glideslope is desirable for 
overshoot characteristics, it also produces an initital tendency for the 
aircraft to pitch up. To avoid this effect, the capture engage logic was 
selected so that the initial tendency of the elevator is to produce a 
negative pitching moment. 

To achieve this property, consider the component of u^> in (107a), 
corresponding to the elevator command. If the glideslope capture mode is 
engaged at the altitude, h , where the elevator command crosses zero, 
it can be seen that the initial elevator command will tend to gradually 
produce a negative pitching moment. 

While the thrust also has a considerable contribution to the pitching 
moment in this aircraft, due to the engine location, the DIALS autothrottle 
and engine dynamics have much longer time constants than those associated 
with the elevator. Thus, changes in the throttle command and thrust are 
of small magnitude over short period of time, so that the initial pitching 
tendency is determined by the faster elevator command. Although the 
initial tendency of the DIALS throttle command is to lower the thrust 
level, which is the desirable action, the glideslope capture criterion 
was selected by considering only the elevator action. After some manip- 
ulation, and neglecting terms of small magnitude, the following glideslope 
capture altitude can be obtained. 
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where Y q is the selected glideslope angle, x g is the estimated ground 
distance to the glide path intercept point (GPIP) , e^ is the estimated 
error in the state x., and z is the commanded inertial angle of attack. 
Note that, as in the localizer capture, the glideslope capture altitude 
is lower when the desired speed is higher. Further note that the criterion 
is applicable to the capture of an arbitrary glideslope by using the 
glideslope angle selected by the pilot as Y in (114) . 


C. LOCALIZER AND GLIDESLOPE TRACK 


As the aircraft gets closer to the localizer or glideslope, the local- 
izer or the glideslope track (or hold) modes are engaged, the order depending 
on the particular approach path. In comparison to the capture modes, the 
main objectives considered in the design of the track modes were to achieve 
quick settling on the localizer and glideslope, insensitivity to wind dis- 
turbances (gusts, shear and steady winds), and accurate tracking of the 
localizer and glideslope. 

Since it is not desirable to start the flare maneuver before the air- 
craft is fully stabilized on the glideslope, a low settling time for the 
track modes is important to enable M close-in captures". The capability 
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to accurately capture and track a shallow or steep glideslope provides a 
flexibility which can be used to avoid vortices when following large air- 
craft, reduce the noise perceived by airport communities, and generally 
increase terminal area operations efficienc}^. 

The localizer track mode is engaged when the aircraft is 30 ft. from 
the runway centerline or 25 seconds after the localizer capture has been 
engaged. Since the main objective of the track mode is to minimize 
excursions away from the runway centerline, the gains feeding back the 
cross-track error estimate to the aileron and rudder are increased gradually 
for a smooth mode transition. Since a type 1 system is desirable in this 
mode, an integrator for the cross-track error is also initiated to avoid 
steady state offsets from the runway centerline. Initially, the integrator 
gains were set equal to a tenth of the gains for the cross-track error, 
i.e., the product of the sampling interval and the cross-track error- 
gain. The gains were then adjusted using simulation results to improve 
the performance under the various conditions expected to occur. 

The glideslope track mode is engaged when the aircraft intercepts the 
glideslope, or 25 seconds after the glideslope capture mode has been engaged. 
The latter criterion is used to ensure that the track mode is engaged even 
if the aircraft does not overshoot and cross the glideslope. To reduce any 
excursions from the selected glideslope, the gains feeding back the deviation 
from the glideslope and the error in the normalized vertical velocity are 
increased gradually using an easy-on. To avoid steady state offsets from 
the glideslope, an integrator on the glideslope deviation is introduced 
when the track mode is initiated. 

During both the capture and track modes, the DIALS filters estimate 
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the components of the wind velocities. These estimates are fed to the 
controller to compute the surface command signals. Thus, to the extent 
that the wind velocities are estimated, the control law performs a gust 
alleviation function; i.e., even if no position or inertial velocity 
errors are present, the wind velocity estimates will produce surface 
commands which tend to counteract the effects of winds, thus reducing 
excursions away from the glideslope. The wind velocity estimates are 
also used to mainatin a constant airspeed. This is obtained by command- 
ing an inertial speed deviation equal to the negative of the wind esti- 
mate in the along-track direction. The commanded airspeed thus remains 
constant at the desired nominal value. 


'£2k 


£ 1 C w£ W £k 


’£2k 


£ i C w£ A w£ W £k 


( 115 ) 

( 116 ) 


D. CRAB AND DECRAB 


When a cross-wind component is present, it is desirable to control 
the yaw angle so that the aircraft is headed into the wind with level wings 
and has no sideslip while remaining on the localizer; i.e., it is desir- 
able to "crab 1 ’ into the wind. If this condition is not specifically 
accommodated, an optimal quadratic regulator will usually both yaw and 
roll into the wind, thus having a non-zero steady state sideslip angle, 
which is not desirable from the pilot’s point of view, while the bank 
angle can be uncomfortable. To obtain the crab condition in DIALS, a 
roll integrator is fed back when the aircraft roll reaches 2°, provided 
that the localizer track mode has been engaged. This ensures a zero 
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steady state roll angle forcing the yaw angle to automatically adjust 
itself to the angle necessary for crab. To obtain crab with a quick and 
smooth transient response, the commanded (aerodynamic) sideslip angle is 
set to zero, i.e.. 



where z_ is the normalized lateral velocity (inertial sideslip) and $ 

j Wo 

is the estimate of the steady wind velocity. As minimizing the sideslip 
is desirable in all the phases of flight, except for decrab, the commanded 
sideslip angle is set to zero at localizer capture initiation and remains 
active until decrab. 

During decrab the aircraft heading is aligned with the runway center- 
line, while the roll angle is used to maintain the aircraft on the localizer. 
The decrab mode is initiated at an estimated altitude of 250 ft. It should 
be noted that this decrab altitude was somewaht arbitrarily selected, and 
can be reduced to initiate the maneuver at a lower altitude. To obtain 
the decrab maneuver, first the commands which produce the crab condition 
are phased out; the roll integrator is phased out and the inertial side- 
slip command z^ is set to zero gradually. Since this corresponds to a 
non-zero commanded aerodynamic sideslip angle, the wind velocity estimates 
fed into the controller produce a decrab tendency. To enhance this ten- 
dency and ensure decrab as a steady state condition, an integrator of the 
heading relative to the runway is introduced to aid the aileron and rudder 
actions necessary for decrab. Initially the integrator gains were set 
equal to a tenth of the heading error gains. These were later adjusted 
using simulation results to improve the performance throughout the flight 
regime. 
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E. FLARE 


During the flare maneuver, the main objectives are to reduce the air- 
craft’s sink rate to an appropriately lower level and to touch down with 
a sufficient margin on the pitch attitude, aligned with the runway, near 
the selected touch-down point. As DIALS is required to flare from steeper 
glideslopes than usual for a B-737, the aircraft has a higher sink rate 
than usual on the glideslope. The approach taken was to generate a flare 
trajectory on line as a function of the glideslope angle, the desired 
touch-down flight path angle and the touch-down point; so that when the 
glideslope is steeper, the flare (initiation) altitude is higher. It 
should be noted that, for a given approach, the flare path is fixed in 
space, and is fed to the automatic control system as the desired or 
commanded trajectory. 

The family of flare paths was developed starting from the desired 
vertical acceleration profile. The following form was selected for this 
purpose 


h"(x ) = 
e 


1+cos 2II(x +AX r ) /P 
e f 

2(H.-X. tan y ) 
f f 'o 


, 2 x +AX^ < P 

e f 


(118) 


otherwise’ 
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AX f = (P - 2 X f )/2 , (120) 


The prime " ’ " denotes differentiation with respect to the distance 

variable x . 

e 


42 



where x is the x-coordinate of the aircraft’s position vector, the 
e f 

touchdown point, the altitude at (zero for touchdown), Y Q and Y ^ 
the selected flight path angles on the glideslope and at touchdown, re- 
spectively. The desirable characteristics of this profile included the 
smoothness of transition (e.g., h u , hence h, has a continuous first 
derivative at flare initiation), the simplicity of the parametric form 
(i.e., flare profiles from various glideslopes are obtained according to 
the value of Y ) > and the case with which h" can be integrated analytically 
to obtain an expression for h’ and h to determine an altitude profile for 
flare. 

It is of interest to note that the flare profile described above was 
selected after some experimentation with an altitude versus x^ profile of 
exponential form. These formulations resulted in undersirable transient 
behaviour at flare initiation, presumably due to the discontinuity in the 
commanded vertical acceleration when starting flare. The cosine type 
acceleration profile provides a smooth transition in this variable, while 
incorporating all the parameters of interest in the flare path such as 
glideslope angle, touchdown point, etc. For a constant ground speed, x^, 
the vertical acceleration profile is seen to be 

h(x ) = x h' (x ) , (121) 

e e e 

h(x ) = X 2 h"(x ) . (122) 

e e e 

The vertical profile h(x g ) , h(x e ), h(x^) is used in the desired path or 
command vectors, z^ and in order to follow the flare path for the se- 
lected glideslope. However, the altitude profile alone does not guarantee 
all of the desirable and critical properties necessary for an acceptable 
touchdown. The pitch attitude at touchdown is of utmost importance. A 
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negative pitch angle at touchdown would result in the nose gear touching 
down before the landing gear and being subjected to high levels of loading. 
Thus, a positive pitch angle safety margin is usually desirable. It is 
also usual to have a decreasing airspeed profile during flare to aid the 
pitch profile during flare, as well as ease the touchdown speed, sink rate 
and rollout. To obtain these characteristics the command variables and 
were set as follows. 
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where Y d ^ is the desired flight path angle at time t^, GEZ5 is an easy-on, 

Z' , is the desired altitude in the earth coordinate axes, L. . is the 
ed k ij 

estimate of L E g(i,j), and is the first column of the identity matrix. 

DELVFR and G are constants which can be set according to the amount of 
^4 

decrease in airspeed and the pitch profile desired, respectively. These 
values were also modified using the profiles obtained in simulation runs. 


F. FILTER DEVELOPMENT 


As seen in (107) , the full state feedback LQG problem requires estimates 
of all the state variables in its implementation. As the states correspond- 
ing to the desired trajectory are known, these variables are not estimated. 
However, the aircraft motion variables and, in particular, the wind velocity 
variables including steady winds, shear winds and gusts are required for 
quick control adjustments to changing wind conditions in the full state 
feedback formulation*. 

The MLS position information and the on-board sensor signals are fed 
to the longitudinal and lateral DIALS filters shown in Fig. 4 and 5. The 
filtered estimates are then used to compute the control surface commands 
as well as to define the desired path. The form of filter used is a dis- 
crete Kalman filter with constant gains. The steady state Kalman filter 
gains were obtained as a starting point, then some gains were modified to 
obtain better performance using a non-linear simulation of the aircraft. 

The use of constant gains for the Kalman filters reduces the considerable 
computation requirements imposed by the error covariance updates. The use 

* 

The author’s recent development of a fast, reliable and convergent algorithm 
for the stochastic output (limited state) feedback problem has made the use 
of optimal dynamic compensation along with complementary filters practical, 
and more desirable in high order problems [»]• 
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of discrete filters along with digital controls is well suited to the dis- 
crete nature of MLS guidance system. The goals of the filter development 
include the accurate estimation of the aircraft’s position, velocity, atti- 
tude and wind velocities using on-board sensors usually available on commer- 
cial aircraft, but without using costly inertial platforms, and angle of 
attack or sideslip sensors which are not currently available on many aircraft. 

The aircraft's position is obtained using the MLS guidance system, 
which provides volumetric coverage in the terminal area. The aircraft re- 
ceives range, azimuth and elevation information at discrete intervals from 
which it can obtain its position with high accuracy even under adverse 
weather conditions. The ground azimuth antenna is located at the runway 
centerline with coverate up to ±60° . The DME antenna which provides the 
range of the aircraft is generally co-located with the azimuth antenna. 

If the DME is located to the side of the runway, a simple transformation 
can be used to obtain the aircraft's position. The elevation antenna is 
located at the glidepath intercept point (GPIP) , but is offset to the side 
of the runway; it provides the aircraft's elevation angle up to 20°. Thus, 
the aircraft has accurate position information in the volume of space with- 
in the limits mentioned above. 

Consider a right handed coordinate frame with its origin at the phase 
center of the azimuth anetnna, the x-axis along runway centerline and pos- 
itive towards the runway, and the z-axis positive vertically upwards. If 
the position of the aircraft in this coordinate system (the MLS coordinate 
frame) is ( X Q > y Q > -Z Q ) » then the MLS signals have the values given by the 
formulas below. 

V ~2 2 2 

x o + y o + z o (135) 
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(136) 



El = tan -1 —• — ■ ■ ° ===- (137) 

V (x o " V + (y o " V 

where and y are the x and y coordinates of the elevation antenna phase 
center in the MLS coordinate frame. 

The onboard sensors used are three body-mounted accelerometers, attitude 
gyros for the pitch, roll and yaw angles, attitude rate gyros, barometric 
altitude and sink rate, airspeed and radar altitude. The filter does not 
require expensive inertial platforms for sensor measurements. All of the 
onboard sensors are usually available on commercial aircraft except for the 
body accelerometers. However, the accelerometers are relatively inexpensive 
instruments and in many of the newer aircraft, a normal accelerometer already 
exists for pitch axis control. It should be noted that the radar altitude 
measurement is used in place of the MLS elevation signal only after the air- 
craft crosses the runway threshold since the aircraft flies out of the MLS 
elevation coverage during the flare maneuver. 

To provide accurate estimates of the wind velocity components in the 
longitudinal, lateral and vertical directions, the filters were formulated 
using the aerodynamic properties of the aircraft. In addition, estimates 
of the bias errors in the three body-mounted accelerometers, barometric 
altitude and sink rate and attitude sensors are also obtained by the filters. 
The lateral filter uses the roll and yaw attitude, yaw rate and roll rate, 
the processed MLS y^ measurement, and the lateral accelerometer. The 
remaining measurements are used for the longitudinal filter. Discrete 
longitudinal and lateral Kalman filters were formulated with steady state 
gains, with some cross-coupling using second order terms and wind estimates. 
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The measurement bias terms were not included in the Kalman filter formula- 
tion and the steady state gain computations, but the required gain matrices 
were obtained using simple approximations. The basic form of the filters 
can be expressed by 


\ = y k - C x S k - c w w k - c b b k > 
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V. INITIAL SIMULATION RESULTS 


As an initial evaluation of the basic approach used in the design of 
DIALS (i.e. , a digital modern control structure tracking a trajectory 
defined on-line), and of the particular concepts used to obtain the various 
characteristics desired in each phase of the final approach and landing, 
the automatic landing system design and a six-degree-of-f reedom non-' 
linear aircraft model were simulated on a digital computer. The results 
of this initial evaluation are given in this section. A more detailed, 
realistic and in-depth evaluation, and further development of the auto- 
matic control law was performed prior to the flight testing of DIALS. The 
latter evaluation used an aircraft dynamics simulation validated by actual 
flight data for the TCV Research Aircraft, and realistic non-linear engine, 
actuator and servomechanism models [10] The results of the detailed 
evaluation will be given together with the DIALS flight test results in a 
separate report. 

The simulation used in the initial evaluation is a non-linear, six- 
degree-of-f reedom, rigid aircraft model. The aerodynamic forces are 
generated using the stability derivatives of the Boeing 737-100 aircraft 
at an airspeed of 61.73 m/sec or 120 knots. The stability derivatives 
remain constant throughout the simulation; however, the dynamic pressure 
and the attitudes used in generating the aerodynamic forces and moments 
vary according to the path and wind conditions. Non-linear kinematic 
equations are used to generate the position of the center of mass and the 
attitude of the aircraft assuming a rigid structure. A detailed description 
of the aircraft equation of motion used in this simulation can be found 
in [20]. 
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The servomechanisms and hydraulic actuators on the aircraft are modeled 
using first order dynamics by introducing an average lag in the control com- 
mands; however, non-linearities such as hysteresis or backlash are not 
included. The engine dynamics is modeled by a first order system whose 
time constant varies depending on whether the thrust is increasing or de- 
creasing. The mechanical link from the aileron to the spoilers is modeled 
as an instantaneous, directly proportional motion as shown in (65). 

Except for the engine dynamics, the overall actuator systems are 
assumed to be linear and contain no rate limits. However, it should be 
noted that all the DIALS control commands contain position or rate limits 
which are below the aircraft's actuator limits and capabilities. Thus, 
the lack of rate limits in the simulation does not introduce important 
errors in the results. However, non-linear actuator effects which are not 
considered part of this evaluation, but will be analyzed later. 

The phases of flight considered in the simulation are the localizer 
and glideslope capture, the localizer and glideslope track, the crab/decrab 
maneuver, and flare until touchdown. Both the usual 3° glideslope and the 
steeper 4.5° glideslopes are simulated and the control performance evaluated 
in weather conditions that include gusts and steady winds. A simulation 
of the final approach that contains shear wind conditions for a steep ap- 
proach (4.5° glidelsope) are also shown. 

To evaluate the initial transient effects of the control law when it 
automatically engages the localizer and glideslope capture modes, the 
simulation is started prior to the engagement of DIALS. Thus, at the begin- 
ning of the simulation, the aircraft is in level flight, approaching the 
runway at the reference airspeed for landing (selected by the pilot) on a 
straight line which intercepts the runway centerline at a specified track 
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angle knd corresponding yaw angle. For the simulation results shown in 
Figures 5 through 8, the aircraft is tracking a 30° intercept path to the 
runway centerline and flying a constant altitude path under the automatic 
control of a 3-D area navigation guidance and control mode. 

This provides a realistic method of analyzing the initial transition 
effects, in comparison to starting the simulation at the instant DIALS 
engages, since it simulates the flight conditions more realistically. In 
particular, note that due to wind conditions and control errors, the air- 
craft does not always have exactly zero bank angle at' localizer capture 
initiation, or exactly zero sink rate at glideslope capture initiation as 
can be seen in Figures 7, 8, and 9. Thus, transition effects can be 
analyzed more realistically. 

It should be noted that similar transition effects are present as 
DIALS automatically engages each hew control mode. Thus, transition from 
the localizer and glideslope capture modes to the corresponding track modes 
followed by transition to the decrab and flare modes can cause undesirable 
transient effects, as in each case, control law changes such as some gain 
changes, introduction of integrators or changes in the commanded path, 
speed or attitude occur when each new mode is engaged. Observation of the 
simulations given in Figures 6 through 9 shows that such mode-to-mode 
transitions occur smoothly and without any undesirable transient effects 
in no wind as well as turbulent conditions. In fact, it would be difficult 
to determine the exact time the mode transitions occur from the plots shown 
The smoothness of the mode-to-mode transitions was obtained using easy-on 
functions for gain changes, and special attention to the selection of the 
desired or commanded trajectory, as described in the section on flare. 
Smooth easy-on functions were not used when initiating or terminating 
integral feedback. 
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Figure 6 shows the simulation of a steep (4.5°) final approach and 
landing at a reference airspeed selection of 64.3 m/sec. or 125 knots in 
a no-wind condition and perfect measurements; i.e., no sensor error, noise 
or bias is present in this run. Thus, this run provides a baseline from 
which the effects of wind disturbances and measurement errors and noise 
can be analyzed and evaluated. These effects can be seen in the simulation 
shown in Figures 7 through 9 which contain both steady winds and gusts as 
well as measurement errors, noise and biases. In particular. Figure 9 
shows the final approach and landing simulation in a wind shear condition. 

The wind gust conditions are simulated according to the Dryden spectra. 
The models used for generating both the lateral and longitudinal wind gust 
velocities are given in Section II. The standard deviation of the wind 
gust velocities for the simulation runs shown here is .61 m/sec. or 2 ft/sec. 
in all three directions, with corresponding levels for moments due to gusts 
determined by the Dryden models. The steady winds are fixed relative to 
the earth coordinate axes. For the simulations shown in Figures 7 and 8, 
the steady wind is a 5.14 m/sec. or 10 knot quartering headwind. The shear 
wind conditions in Figure 9 simulate a steady headwind of 10 knots at 
altitudes above 200 feet. Between 200 feet and 100 feet, the headwind 

decreases at a rate of 4 knots per 100 feet. Below 100 feet, the headwind 

decreases at a rate of 8 knots per 100 feet. Thus, a 10 knot headwind at 

200 feet changes to 2 knot tailwind at touchdown. 

Depending on wind conditions, the localizer capture mode engages 
between 10 to 15 seconds after the start of the simulation. The exact 
time of the capture initiation is determined by the captive criterion 
given in (113) . The baseline run in Figure 6 shows an initial inclination 
to roll away from the runway (as desired), reaching a maximum of about 8° 
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bank angle before settling at a level wings conditions. The cross-track 
error or distance from runway centerline has an overshoot of approximately 
40 ft., and settles in approximately 50 sec. after capture initiation. 

The aircraft heading is similarly aligned with runway in about 55 sec. 
after the capture. The effects of the quartering headwind (corresponding 
to 7 knot cross-wind) , as well as the effects of sensor noise and errors 
are seen to produce small deviations from the baseline response such as 
to increase the cross-track error overshoot to a maximum of about 75 ft. 
and to produce an initial tendency to roll towards the runway before re- 
versing the roll direction due to initial errors in the estimates, and 
approximations used in developing the capture criterion.* However, the 
settling time is seen to be largely unaffected by wind conditions and 
remains in the vicinity of 50 seconds. In comparison to current ILS systems 
using inertial platforms and classical design structures, the overshoot 
performance shown here is of similar magnitude; however, the settling time 
required for DIALS is less than half the time required for these systems 
as can be seen in Figure 10 showing the response of a typical autoland 
system. 

The glideslope capture mode is initiated between 22 to 27 seconds 
after the start of the simulation and vary shortly after the localizer 
capture mode is engaged. The variation in the capture initiation time 
is due to differing wind conditions and to different glideslopes. As can 
be seen from the pitch angle and sink rate behaviour prior to capture, 
the aircraft is not generally settled and trimmed on a constant altitude 
path when the glideslope capture mode engages. However, the capture mode 

k 

Both the initial rolling tendency as well as the overshoots were consider- 
ably improved in futher development and will be shown in the flight test 
results . 
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produces an initial pitch down action in each case. Similarly, the sink 
rate is quickly reduced from its initial level to that required to follow 
the selected glideslope at the selected airspeed. It should be noted that 
there is practically no overshoot in the sink rate. The glideslope devi- 
ation similarly is reduced quickly. It should be noted that the aircraft 
is within 10 ft. of the glideslope within 15 seconds after the capture 
mode is engaged, and does not overshoot the glideslope. Finally, the 
capture of the steep (4.5°) glideslope is performed without any apparent 
difficulty. The effects of the wind conditions and sensor errors seem to 
produce minor deviations from the baseline run. 

The localizer and glideslope track modes can be seen to provide the 
necessary control action to remain on the localizer and glideslope without 
any appreciable offset in different wind conditions. The type-1 property 
provided by the integral feedback appears satisfactory. It should also be 
noted that the deviations from the localizer and glideslope due to wind 
gusts tend to be small and are quickly corrected due to the gust alleviation 
or feedback of wind velocity estimates in the control law. 

In all the simulation runs, the localizer capture mode is engaged 
first, followed by the glideslope capture within a few seconds. Thus, 
both capture maneuvers occur simultaneously, and settle on the desired 
path within 50 seconds. This simultaneous action and quick settling 
enables the capture maneuvers to start closer to the runway in the per- 
formance of close-in captures. The ability to use a shallow or steep 
glideslope according to termainal area conditions provides a further 
flexibility that can improve terminal area operations. 

When a cross-wind is present, DIALS automatically crabs the aircraft 
into the wind while maintaining level wings, as desired. When the decrab 
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altitude is reached, DIALS performs a sideslip maneuver to align the air- 
craft heading with the runway. As can be seen in Figures 7 and 8, this is 
achieved by rolling the aircraft into the wind while reducing the heading 
relative to the runway. 

The last control mode initiated is the flare maneuver. For the 3 and 
4.5 degree glideslopes, the maneuver is initiated respectively at 31.2 
meters (102.2 feet) and 40.8 meters (133.7 feet) of altitude. In all the 
simulation runs, the desired touchdown point is 396.2 meters (1300 feet) 
from the GPIP and the desired touchdown vertical velocity is .67 meters/ 
sec (-2.2 ft/sec) at 64.3 meters/ sec (125 knots) reference airspeed. 

For the -4.5° glideslope no wind no noise case the touchdown point and 
touchdown vertical velocity were respectively 381.9 meters (1253 feet) 
and -.88 meters/sec (-2.89 ft/sec). 

For the cases of sensor noises and winds, the touchdown points and 
vertical velocities for the 4.5° glideslope with no shear, 3° glideslope with 
no shear, and 4.5° glideslope with shear are, respectively, 353.3 meters 
(1159 feet) and .79 meters/sec (-2.58 ft/sec), 377.7 meters (1239 feet) 
and .59 meters/sec (-1.93 ft/sec), and 436.8 meters (1433 feet) and .98 
meters/sec (-3.29 ft/sec). For these cases the touchdown point is 
within 43.0 meters (141 feet) of the desired point. The highest touchdown 
sink rate, occuring for the shear wind case, exceeded the commanded sink 
rate by .30 m/sec (1 ft/sec). It should also be noted that in all the 
cases the touchdown pitch attitude is sufficiently positive to prevent 
premature nose wheel touchdown. The position offsets from runway center- 
line for these simulation runs were, respectively, .39 meters (1.27 feet), 
1.22 meters (4.00 feet), .18 meters (.59 feet), and -.41 meters (-1.34 
feet) . 
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It is of interest to note that the wind shear is encountered while 
tracking a steep (4.5°) glideslope. As the wind velocity variation is 
simulated according to aircraft altitude, the time rate of change of the 
wind velocity is higher when flying a steeper glideslope. Thus, the 
effect of the wind shear on a 4.5° glideslope is 50% higher than on a 
3° glideslope. Nevertheless, Figure 9 shows that the automatic landing 
system handles this wind shear scenario without any difficulty, and 
tracks the steep glideslope and flare path without introducing any extra 
altitude error. 
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VI. CONCLUSIONS 


The design and development of a digital automatic landing system for 
a small commercial jet transport using modern digital control techniques 
has been considered. The system uses MLS position information, body- 
mounted accelerometers and on-board sensors usually available on commer- 
cial jets, but does not require costly inertial platforms. The phases 
of final approach to landing considered in the design were: 


1) 

Localizer 

capture 

2) 

Steep Gideslope capture 

3) 

Localizer 

track 

4) 

Steep Glideslope track 

5) 

Decrab 


6) 

Flare. 



The simulation results presented indicate that the digital integrated 
automatic landing system can provide improved performance, without the 
use of inertial platforms, by providing 1) added flexibility in capturing 
and tracking steep glideslopes, 2) reducing the settling time to enable 
close-in captures, 3) reducing the sampling rate to 10 Hz to reduce compu- 
tational requirements, 4) estimating wind velocities to provide better 
path control in wind shear conditions. The capability of flying a select- 
able steep glideslope, low overshoots and short settling times indicate a 
positive impact on airspace utilization, terminal area capacity in adverse 
weather conditions, perceived noise in airport communities, and avoidance 
of wake vortices. 

The results have shown that the application of direct-digital-design 
and modern control theory methodologies can be used to engineer an 
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automatic landing system that performs the various lateral and longitudinal 
maneuvers with minimal gain changes using a fixed control law structure. To 
the author's knowledge, the digital integrated automatic landing system 
(DIALS) described in this report is the first automatic landing system 
designed using modern digital control methodologies whose performance has 
been demonstrated in successful flight tests. 
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FIGURE 1 DIALS FLIGHT PATH GEOMETRY 















FIGURE 4 DIALS FUNCTIONAL BLOCK DIAGRAM 



FIGURE 5 


BLOCK DIAGRAM OF FEEDBACK LOOP 
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FITURE 7 4.5° GLIDESLOPE, STEADY WINDS, GUSTS, SENSOR NOISES 

AND BIASES 
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FIGURE 10 RESPONSE OE AN INERT IALLY AUGMENTED ILS AUTOLAND 
SYSTEM 3° GLIDESLOPE, STEADY WINDS, GUSTS, SENSOR 
NOISES AND BIASES 
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